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CONSTRUCTIVE  FUNCTION  THEORETIC  METHODS  FOR 


HIGHER  ORDER  PSEUDOPARABOLIC  EQUATIONS* 

R.P.  Gilbert  and  G.C.  Hsiao 

0.  INTRODUCTION 

In  this  work  we  will  develop  a constructive  method  for  solving 
pseudoparabolic  equations  of  order  2n  in  the  plane.  More  precisely, 
we  investigate  equations  of  the  form 


(0.1) 

«£lu]  :=  Mlut)  + 

L[u  ] , 

where 

M and  L are  the  respective  elliptic  operators 

n 

(0.2) 

M[u]  :=  Anu  + 

l M <An'ku) 

k 

= 1 

p+qi  k 

MkUl  :=  1 

p,q=o 

pq  . ' 3P+q» 

aPMx.y  — , 

k 3xp3yq 

and 

(0.3) 

m 

L[u]  :=  Amu  + 

L^(Ara  ku)  , m < n 

k=l 

p+q<k 

LkUl  - f b£q 

p , q-o 

(x,y)  — - — - . 

3xp3yq 

The  coefficients  of  , L^  are  taken,  furthermore,  to  be  analytic 

functions  of  x and  y for  (x,y)  € D C.  C^. 

Recently,  the  integral  operator  methods  of  BERGMAN  ( 2 ] and 
VEKUA  [22),  which  have  been  very  successful  for  developing  repre- 
sentations for  solutions  of  elliptic  equations  in  the  plane,  have 
been  extended  by  COLTON  to  treat  the  cases  of  parabolic  [ 9 ] and 
second  order  pseudoparabolic  equations  [10],  [111  in  the  plane. 

BROWN,  GILBERT  and  HSIAO  [ 8 ] and  BROWN  and  GILBERT  [ 7 ) developed 
analogous  techniques  for  fourth-order  pseudoparabolic  equations  using 
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respectively  the  methods  of  VEKUA  and  BERGMAN  to  treat  elliptic 
operators.  BROWN  ( 6 ] , on  the  other  hand,  completed  the  study  of 
fourth-order,  analytic,  parabolic  equations  in  two  space  variables. 

Investigations  concerning  integral  operators  which  generate 
solutions  to  parabolic  and  pseudoparabolic  equations  in  three  and 
four  space  dimensions  have  been  made  by  RUNDELL  [12],  RUNDELL  and 
STECHER  [13],  and  by  BHATNAGAR  and  GILBERT  [ 3 ],M  4 ],  [ 5 ]. 

Pseudoparabolic  equations  arise  in  a variety  of  physical  problems, 
such  as  the  velocity  of  a non-steady  flow  of  a viscous  fluid  [21], 
the  theory  of  seepage  of  homogeneous  fluids  through  fissured  rock  [1], 
hydrostatic  excess  pressure  during  the  consolidation  of  clay  [19], 
and  the  stability  of  liquid  filled  shells  [18],  [23],  [24].  GILBERT 
and  ROACH  are  presently  investigating  the  last  mentioned  problem  as 
an  application  of  some  of  the  ideas  presented  in  this  current  work. 

I.  THE  FUNDAMENTAL  SOLUTION 

If  the  coefficients  of  as  analytic  functions  of  x,y 

have  an  analytic  extension  to  (z,z*)  £ D x D*  , D*  :=  {z:  z 6 D)  , 
where  z = x + iy,  z*  = x - iy  , then  M[u]  has  a representation 


(1.1) 

M[u] 

n 

:=  £ A . 

k , j=o  k3 

(«,**)  3k+jU  , 

szkaz*3 

U ( z , z* ) 

* # 

/ Z + Z Z- z \ 

u(--'  -- r) 

with 

Ann  5 

1.  We  assume 

also  that  the  operator 

L )ias  a complex 

form 

(1.2) 

L[U]  : 

n- 1 

“ l B (z,z#) 

k , 3=0  k;) 

?k+3u 

3zk3z^ 

The  adjoint 

operator  to 

X is  given  by 

X*  =/f, 

(1.3) 

/flu) 

:=  M*[Ut]  - L* [ U 

] , 

with 

(1.4) 

M*  [ U 1 

:=  [ <-l)k+3 

k , j=o 

3k^ 

(AkjU)  , 

3zSz#J 

(1.5) 

L*[U] 

n- 1 . 

:=  l <-l)k  + :) 

k , 3=0 

3k  + ^ 

(hmu)  . 

3Zk,,Z*] 
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A function  S of  the  form 

(1-6)  S (x,y , t ; f, , n . T ) : =*  A(x,ylt;(,n<T).  In  - + B(x,y,t;C»n,T) 

2 2 W 

where  r = [(x-f.)  t (y-rj)  )T  will  be  called  a fundamental  solut  ion 
of  (1.1)  if  it  satisfies  the  following  conditions 

(c-1)  As  a function  of  (x,y,t),  S is  a solution  of  the  adjoint 

equation  .Al  (S ] = 0 and  is  an  analytic  function  of  its  argument 
except  at  r = 0,  where  An  has  a logarithmic  singularity. 

(c-2)  At  the  parameter  point  x = t,  y = n we  have 
3P+qA 

= 0 for  p+q  <.  2n-3  , and  An  ',\  = 4n  | 

3xp3yq  L 

(c-3)  A and  B are  analytic  functions  of  (x,y,t)  at  r = 0 

and  vanish  at  t = T . 


Remark : The  above  implies  that  A may  be  written  as 

2 n — 2 n a 

r A (x,  y , t ; £ , n , f ) witli  A regular  at  the  parameter  point.  The 
above  notation  is  computationally  easy  to  work  with. 

We  intend  to  show  that  it  is  possible  to  develop  the  coefficients 
A and  B as  analytic  functions  with  the  expansions 


(1.7) 


(1.8) 


A(z,z*,t;C,C*,7) 


B(z,z*,t;C,C*,T) 


l A. (z,z*;C,t*) 
j = l -1 


l B. (Z,Z*;C,C»> 
j = l 3 


(t-f ) 
3 ! 


( t-T ) j 


3! 


Furthermore,  we  shall  identify  A^  (z , z * ; C , C* ) as  the  Kiemann 
function  corresponding  to  the  operator  tj  • The  other  coefficients 
will  be  seen  to  satisfy  homogeneous  Goursat  conditions  on  z = C 
and  z*  = !.*  . 

Insert  ing  N(  1 . 6 ) into  Al  ( U)  :=  MMU^J  - L*(U1  T;  0 we  obtain, 
after  some  manipulation  of  terms 

(1.8)  /t[S|  H K (Al  in  i + #(B)  + I + I * = 0 

/V  V *•  ^ I l * 1 


* 
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with 


a. io) 


n- 1 

+ y 

[>  ) 

Vi)k+n 

3k+n-p 

k=p 

\k-p/  | 

3zk~p3z*n 

n-1 

+ l ( 

:-i)k+j 

gk+ j -p 

-r  (A,  .A  • 

j=o 

3zk~P3z* 

] k3  t 

■ (A  .A  J 
n-p,  *]  "3  t 


A,) 


and 


I* 

n 


j , ■ r 


I'n-D!  jo(-DkS 


3zk 


r^(:.q)j  <-D 

q - 1 \ ^ / k=o 


k+n  3k+n-q  . 

3zk3z.n-q  kn 


(1.11) 


(-Dn+j  -i 


n + ]-q 


(A  . A,  ) 
3zn3z.l-<l  n3  t 

(A.  .A. 


+ ny1(-i)k^j  3**3"?  q (A.  .A.  - B.  ,A)1  ).<--A)q..(S:ALL  . 

3zk3z*3_q  kj  t k]  jj'  (z*-C*)cl 

Because  of  the  raultivaluedness  of  the  logarithmic  singularity  it  is 
necessary  to  set  A [A]  = 0 . To  cancel  the  poles  at  z = C and 
z*  = c*  we  ask  that  the  coefficients  of  (z-^)p  , (z*~C*)p  , 

(p=l<...»n)  vanish.  These  latter  conditions  provide  us  with  so- 
called  Goursat  data  for  the  A^ , B ^ coefficients  in  the  repre- 
sentations (1.9).  Setting  (1.7)  into  (1.10)  and  (1.11)  we  obtain 
the  following  conditions  which  A^  (z  , z*  ; t; , 5*  ) must  satisfy 


(1.12) 


Hi)  l -4^7  kj  A ( C , z*  ; C , t*,l  - 

k=p\P  / j=o  3zK  PSz*:>  L k3  1 J 


(P=l,2, . . . ,n)  , 


and 
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(1.13) 


X [I)  l <-l>k+3  3k  3 j-q  [Akj(z,C*)  A (z,;*;t,C*)l  = 0 
3=q  ' 'l  / k=o  Sz^Sz*3  q L J J 


(q=l,2, . . . ,n) , z*  = c*  . 

We  recall  that  condition  ( c — 2 ) implies  that 


(1.14) 


3z^3z*^  ^ 


A. (C,C*;C,C*)  = 0 , for  j=l,2,3,...  , 


(p, q=0 , 1 , 2 , . . . , 2n-3  , with  p+q  £ 2n-3)  . 


In  order  to  show  that  A^  is  actually  the  Riemann  function,  it 
is  sufficient  for  us  to  show  that  the  above  conditions  are  equivalent 
to  the  characteristic  conditions  which  uniquely  determine  it.  We 
recall  from  VEKUA  [22),  Chapter  V,  the  following  conditions  imposed 
on  the  Riemann  function  R (z , z* ; 5 , 5* ) and  for  convenience  we  label 
these  conditions  using  his  equation  numbers  : 

(37.28)  M* [R]  = 0 


(37.29) 


(37.30) 


3 ^ R 3 ^ R 

(z,z*;C,C*)  = 0 , (Z,z*;c,;*) 

3z  . . 

z=;  ^ z*=t* 

(k=0 ,1,2, . . . , n-2 ) 

1r 

< z , z * ; C , C * ) = X(z*,c  *,;)  , 

3z 

z=C 

^ 1 D 

- — (z.z*;C.t*)  I = X*(z,C,5*)  . 


Here  X,  and  X*  are  solutions  respectively  of  the  ordinary  dif- 
ferential equations 


n , ,k  r 

* (-1}  A"kU'2 
k=o  3z* 


l (-1)’ 


Ank(C,z*)  X = 
Akn(2,C*)  X* 


= 0 , 


Furthermore,  X and  X*  are  seen  to  satisfy  the  intial  conditions 


(37.32) 


( z * , t * , r. ) 


( z * . c * . ‘r. ) 


( k-0  , l,  . . .n-2) 


3n-lX* 


It  u'r-'r-*)  I - 0 - 


( k=o , 1 , . . .n-2) 

We  note  that  ( z , z*  ; c , c.  *)  automatically  satisfies  ( 37.28), 
(37.29),  (37.32)  by  virtue  of  our  conditions  (c-1)  and  (c-2) . This 

suggests  that  we  check  to  see  if  the  remaining  conditions  (37.30), 
(37.31)  are  compatible  with  ours.  Using  (37.29)  in  combination  with 
condition  (1.12)  we  obtain  for  p > 1 

l (-Dn^  f)  - ^ 


(1.15) 


j=o 

Dz 

case 

where  p 

= 1, 

n 

1 <- 

, • n- 1 
un+3  l 

('V 

j=o 

l=o 

pn:'2’2*'  Ai(2'2*^'r-*>  j 


(1.16) 


n+i  "T-1  / n-l\  3 ) 3 

11  i l ^ 'j 

j0(-1,n+i  \w^  h^r  <z'2*“'‘*H  ] - 0 • 


Identifying  temporarily  A ^ with  R and  identifying  the  (n-l)s 
derivative  with  respect  to  z as  X as  given  in  (37.30),  (1.16) 
becomes  the  first  of  equations  (37.31).  Repeating  this  analysis 
with  (1.13)  and  identifying  X*  as  the  (n-l)st  derivative  with 
respect  to  z*  , (37.30)  yields  likewise  the  second  equation  of 

(37.31) .  This  exhausts  our  conditions  and  leaves  us  free  to  impose 
the  additional . initial  conditions  on  X and  X*  prescribed  by 

(37.32) .  We  conclude  that  we  are  permitted  to  identify  Aj  R . 

Lemma : The  f i rst  coefficient  A<  (ZrS*f  C4C*)  for  A i_n  the 

representation  (1.7)  may  be  taken  as  the  Rlemann  funct ion  associat< 
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with  M [ u ] = 0 . The  coef  f icients  A(z,z*;c,c*),  p i 2 , are 
determined  recurs  i ve ly  by  the  nonhomoqeneous  equa  t ions 
(1.17)  M*  [Apl  = L*  (Ap_1)  , 

and  the  homogeneous  Goursat  data  , 


(1.18) 


A (z,z*;C,c*) 
3 z H 


= 0 


z=t 


A 


P 


(z , 


z* 


1=0,1, • • • , n-l 


N 

1=0,1, n-l  . 


Proof : Since  S = A log  i + B satisfies  the  adjoint  equation,  it 

follows  that^llA)  = 0.  This  in  turn  implies  M* [A^ ] = 0 , and 
M*(Aj  + ^]  = L*  [Aj ) , ( j = 0 , 1 , . . . ) . A moment’s  reflection  concerning 

our  condition  ( c— 1 ) indicates  that  it  can  be  satisfied  using  the 
above  conditions  (1.18),  when  A^  is  taken  to  be  the  Riemann 
function  for  M|U]  = 0 . 


II.  DETERMINATION  OF  THE  COEFFICIENTS  A (z,z* ;;,;*) 

The  series  representations  (1.9)  for  the  coefficients  A and 
D of  the  singular  solution 

S(z,z*,t;c,t*,t)  :=  A (z , z * , t ; ; , t* , t ) In  i + B (z , z* , t ; C , , i) 
suggest  that  we  try  to  determine  Ap  and  B successively.  To 
this  end  we  develop  a Green's  formula  based  on  the  formal  identity 

(2.1)  V M[U)  - U M*  [ V ] = ~(z,z*)  + |^(z,z*)  , 

where  M and  M*  are  given  by  (1.4)  and  (1.7  and 


(2.2) 


and 


P*  <Z,Z*) 


n 

y 

k-l 

y 

1 

3-0 

L 

p=o 

k = l 

n 

j-1 

y 

V 

L 

L 

k -o 

Xi 

II 

0 

,-nP  1! 


(VA„J 


3k+j-p-lu 


3z* 


kj  3zk-p-ls,,l  ' 


k+q 


3z 


(VA.  ■ ) 


kj 


j-1 


9z  * 


i1''5  'u 

, ..  i-q-i 
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If  ys  is  meant  as  the  M operator  with  s,  s*  replacing  z,  z*  , 
then,  setting  U(s,s*)  :=  Ax (z , z* ; s , s* ) and  noting  that  A,  is  also 

a solution  of  Ms [ U ] = 0 , the  identity  (2.1)  yields 

(2.4)  - U ( s , s * ) M*[V]  = ||  ( s , s * ) + |~  (s,s*)  . 

Lemma  2 : The  coef  f icients  A ^ (z  , z * may  be  f o rma  1 1 y deter- 

mined by  the  recurs i ve  scheme 

Ap+1  (z,  z*  ; t; , c*  ) = j ds*  ( ds  (z  , z*  ; s , s*  ) L*  [A.^  ( s , s*  ;;,;*)  ) , 

(2.5)  C*  C p= 1 , 2 

Proof : Integrating  (2.4)  yields 

- [ ds*  f ds  U (s , s* ) M*[V]  = f2  P (z , s* ) ds*  4 [Zp*(s,z*)ds 


- Q(z,Z*, -£,(_*)  / 


r z * / 2 

(z,z*;e,c*)  = P P(C,s*)ds*  - | *"  P*  (s , 4 


Using  the  conditions  (37.29)  and  recalling  U :=  A^  , we  can 

simplify  the  integral  of  P above  as 

’ z * n r 2 * n-l  + j 

(2.6)  P (z  , s* ) ds*  = Y V (z , s * ) A ,(z,s*)  2 -U(z,s*)  d . 

3=0  ^ OzH-^S*3 


According  to  VEKUA  (22)  we  may  identify  the  function 


(37.10) 


g (s  * ; z , z * ) :=  -p  A.  (z,z*;z,s*) 


which,  furthermore,  satisfies  an  ordinary  differential  equation 


(37.9) 


,n  n-1  ,m 

+ l A (z,c) 


Using  (37.9),  the  right-hand  side  of  (2.6)  is  seen  to  vanish 
identically. 
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We  next  turn  our  attention  to  the  P*-integral , and  note  that 


(37.29)  implies  that 


3z 


--  -.-  ( s . z » ) = 0 , except  for  q = 0,  j = n. 


.i-q-i 


Hence 


(2.7) 


P*  (s , z* ) ds  = t (-1) 
k=o 


z ,k 


3s 


k [V(s'z*>  Akn(s’z*>] 


3n_1  0 (s  , z * ) 


B z 


kn-l 


ds 


The  associated  functions 


, n-1 


(37.10) 


g*  (s;z,z*)  :=  — ^ 

3s* 


A.  (z , z*  ; s , s * ) 


are  known  (24)  to  satisfy  the  ordinary  differential  equation 

n-1  jm_ 

"mn  ' 


(37.9) 


ds  m=o 

and  the  initial  conditions 


,n  * n-i 

i-3_  + i a (s,z*>  3_a_  = o , 

n l*  mn  HI 

ds 


ls;z,z*)\  = 0 , U=0,1,  . . . ,n-2) 


3s 


(37.11) 


iH-1  , 

— f--  (s;z,z*) 

Ss"'1 


= 1 


Consequently,  after  some  regrouping,  it  may  be  seen  that 


z n ,z 

P*(s,z*)ds=  V 
k=l 


3 

3 s 


V (s , z*  ) Akn(s,z 


•’] 


(2.8) 


3 q*  (s  ; z , z*) 
3s* 


(-1 ) 2n— 1 V(z,z*)  A (z,z*)  ~r-T  g*(z;z,z*)  + Q* (5 ;z, z*) 
nn  n j. 


ds 
,n- 1 


= - V ( z , z* ) + Q*  u;z,z*) 


where 


(2.9) 


n k-1 


q * ( c ; z , z * ) ==  - l l (-U*+k 

k=l  l=o  3sK 

3*g*  (r,  ; z,  z») 
34l 


k-l-t 


V(c,z*)  Akn(r,,z*) 
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Putting  together  the  above  terms  we  note  that  the  following  repre- 
sentation for  V(z,z*)  has  been  obtained 


(2.10) 


V (z , z*)  = Q*  (c ! 2 , z* ) - Q(z,z*;C,C*) 


+ ds*  ds  A^  (z  , z*  ; s , s* ) M*[V(s,s*)] 


Setting  V(z,z*)  :=  Ap+ ^ (z , z* ; t ) and  recalling  the  initial  con- 

ditions (1.18),  we  conclude  Q*(i;;z,z*)  5 0 . The  expression  for 
Q (z , z* ; ^ , £ * ) may  also  be  seen  to  vanish  identically  by  virtue  of 


the  identities 


z*  n k-1  „ fZ 

P(C,s»)ds*  = l l (-1) 

j=o  i=o  J 

^ k=l  g 


A,  (C,s*;t 


e 1 'Vi 


,k+j-£-l 


az^-l_13s*3 


,t»)  Aki(C,s») 


U(f,  ,s*)  0 , 


r z n j-1  . 

P* (s , c* ) ds  = l £ (-1 ) 
k=o  q=o 

C j = l 


J 8z*qask  \ 

\ aj-q-i 

i,C*)  r 

/ Dz*3-q-i 


( s , r,  * ; f, , r,  * ) 


_I  u (s , ;* ) ds  E 0 


Our  Lemma  is  proved  at  this  point  by  recognizing  that 
M‘tAp+1>  = LMAp)  . 

Lemma 3 : The  coefficients  A^+  ^ ( z , z * ; r,  , {,*)  may  be  computed 

recursively  by 


(2.11)  A +.(z,z *;^,C*)  = ds*  ds  A (z,z*;s,s*)  L* I A,  ( s , s* ; t, C* ) ) / 
P 'c*  J ^ P 

p= 1 , 2 , . . . . Furthermore , the  series  representation  for 

2 

A (z , z*  , t ; 5,  C*  i T ) converges  un i f ormly  i n the  doma i n (G  x G*  x T) 
where  G i_s  the  domain  of  regularity  of  the  coef f icients  i n the 
x-y  coordinates  and  T i_s  a disk  in  the  complex  t-plane. 


Proof:  Recalling  the  form  of  L and  L*  , we  note  that  we  may 


writ  e 
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Ax  (z,  z*  ; s,  s*)  L*  [Ap  (s,s*  ; ;*)  ) - Ap  (s , s*  ; ; , t,*  ) L(A^(z,z*;s,s*)] 


m~X  ^lL<'i)£i?r[Ap(s's*;t,c*)  B^(s's‘)]^^^Ai(z'z*;s'S‘) 


k+ j-t-l 


■ l (-l)k-^-V  l ( - 1 ) — ^r— k-  A 0(s,s*;c,C*)B  . (s,s*)  i 
k=o  3s*]m=o  as*m3skLP  3 Jss*3'"1'1 

j=l  1 


j~A, (z,Z*;s,S*) 


Integration  gives  the  following  identity 


ds*  | ds  A ^ ( z , z * ; s , s * ) [A  (s,s* ; 5, ) J 


t*  \ 


ds*  f ds  A^(s,s*;  ;,(;*)  L [A^  (z  , z*  ; s , s*  ) ] 

t*  X 

+ H + H*  , 


n-1  k-1  „ fZ* 

H ■■=-  l I (-D  ds*/-^ 

j=o  i=o  3s 

k=l  V*  1 

(2.12) 


k+j-t-1 


*\_\rAD(s's*;c’c*)  Bki(s's*>i  ak-,~-i~ 
Ss^L  p k]  J 3sK  1 1 


. A(z,z*;s,s*) 


n-l  . j-1  rz  I ,m+k  r _ 

H*  :=  l (-D  l <-l)m  ds/— — ^ A <s,s*;t.C*)  B (s.s*) 

k=o  m=o  3s*  3s  [ r J 

j = l « 


(2.13) 


A.  (z , z* ; s , s* 


Using  the  conditions  (1.18)  for  A j and  (37.29)  for  Aj  , the 
terms  H and  H*  are  seen  to  identically  vanish.  This  establishes 
the  recursive  definition  (2.11).  To  show  that  the  series  for  A 
converges  we  note  that  according  to  VEKUA  (22],  [>.  186  A . ( z , z * ; t , r,  * ) 


is  dominated  by 


. ,n-l  1 . . ,n-l 

I z-tl  lz*~  * I 

[ (n-l)  ! )2 


in  (G  *G*) 


G*]2  . 
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Consequently,  using  (2.11)  it  follows  by  induction  that 


(2.14)  | 


c l*-a^:V-rrp-2  „ LIA 

^ I f nin.  O \ I 1 1 


[ (n  + p-2)  ! ]* 


V.  ) II  :=  sup  |L[A,]  I 
(G  x G*)  2 

C :=  sup  |A. | . 

(G  x G*)2 


From  (2.14)  it  is  clear  that  the  series  for  A converges  in  the 
stated  domain. 


III.  DETERMINATION  OF  B(z,z*,t;C,C*,t) 

From  (1.10)  we  obtain  as  the  partial  differential  equation  for  B, 
(3.D  Aim  = - in  - in*  , 

where  In  and  I * are  given  by  (1.11)  and  (1.12)  respectively. 

Putting  the  series  expansion  (1.8)  for  B into  the  left-hand  side 
of  (3.1)  yields  then  a recursive  scheme  for  the  coefficients  , 

namely 

(3.2,  MMB.1  = f0(z,z*;c.C*>  =--i  l ^<^-1  l (M 
1 2 p=l  (z-;)P  | k=p  \ Pl 

■ 1 >,l 

■ k-p.  .1  ki  1 / 

j=o  3z  ' 3z*  J 


1 \ P\  i (-Dk+3 

2 q=l  (z‘-c*)q  |j=q  \CW  k=o  3zk3z*J 


(A.  . A,  ) \ , 

3zk3z*J'q  1 / 


(3.3)  M*[Bt+1]=f(|Z,z';C,t*):=I*[B,)-^^  I(-l)j-^. 


(z-On  j = o 3z*l  "3  t+1 


i n7-'  (-l)P(p-l)i  3 (M  ? . nk+j  3k+3-P  . . 


n- 1 / k \ n- 1 ...  -,k+ j-p 

M l <-nk+]  -nti— t <BkiAt,  ♦ 

k p V P / j o 3zk  1 3z*]  kl  1 
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1 (n-1)  ! 

2 


(z*-;*)  k=o 

- J "l  i 

* q=l  (z*-c*)q 


n k 3k 

l «“1>  A 


3z 

n 

l 

j=q 


kn  * + !' 


n W4--4  ak+j"q 

a.-1' 5 


0.0, 


- "j1  PI  V,-,."*! 


j =q  \ Q / k=o 
*=1,2,...  . 


[ (-D’ 


k (BkiV 

32  3z,]  q k;!  *• 


If  we  specify  that  the  satisfy  the  homogeneous  Goursat  data 


(3.4) 


3kB£(z,z*  ;;,;;*) 
3kBt (z,z*; c, C*) 

r~ k 

3z 


z*=4* 


z=S 


= 0 


= 0 


k~0 , 1 , 2 , . . . , n-1 

£=1,2,...  , 


then 


B£+l  (z,z*;4,C*)  = 


ds* 


ds  f (s,S*;C,t*)  Ai(z,z*;s,s* 


£=0,1,...  . <> 

The  majoration  for  the  > while  technically  involved,  proceeds 

by  the  usual  methods  [10],  [11],  (8  ). 

We  summarize  the  discussion  of  Sections  II  and  III  in  the  fol- 
lowing theorem. 

Theorem  1:  Assume  that  the  coefficients  A.,  , B.,p  are  analytic 

]K  - 

functions  of  two  complex  variables  z,  z*  in  the  bicylinder  D*D‘. 
Then  A(z,z*,t;C,t*,t)  and  B(z,z*,t;C,4*,T)  are  analytic  functions 
o f their  six  i ndependent  variables  for  all  ( complex)  t , t , and 
2,;  f D , z*,(.*f  0*  . Moreover , both  A and  B can  be  repre- 

sented by  a tin  i tormly  convergent  series  expansion  in  the  form  of 
(1.7)  and  (1.8)  respect  I vely . 


IV.  BOUNDARY  VALUE  PROBLEMS 


Wo  now  proceed  to  use  the  fundamental  solution  il.B)  t develop  re- 


presentations for  boundary  value  problems.  The  'lot  > us  to  ! on  oity 


value  problems  we  will  seek  in  the  class 


(4.1)  -f f :=  { u ( X , y , t ) : uf  C2n_"  (D  * T)  ; U(.  f C/n  (D  « T)  D U“'n‘  ‘ (D  « T)  ) , 


where  T :=  (t  : 0 < t < t } and  where  t is  a fixed  constant. 

o o 


We  begin  by  considering  the  identity 


't«£(uJ  - ut/((v]  H (vt  M|utJ  - ut  M*[vt)} 


(v  L(utl  - ut  L*  I v]  ) + y-  {v  Udl) 


If  we  replace  v by  the  fundamental  solution 


S(x,y,t;£,n»i)  = A log  — + B , 


then  (4.2)  leads  to  an  integral  representation  for  solutions  which 


satisfy  the  homogeneous  initial  data  'U(x,y,o)  = 0 in  D,  namely 


(4.3)  u(£,n»T)  = 

2n[  (n-1)  ! r 


| dt  [ H [ u ( x , y , t ) , S (x  ,y  , t ; C , n < t ) ] 

> n on 


3k+jut  3k+jS 

Here  11  is  a bilinear  form  in  — r — -r  , — r , (k+j  s.  2n-l)  as 

3x  3y*  3x  Dy4 


3k+ju  g k + j $ 

well  as  — n , — r — -r,  (k+j  <_  2n-3)  . The  introduction  of  a 

3x  3 y ^ 3x  3 y 4 


special,  Green-type,  fundamental  solution,  namely  one  which  satis- 


fies on  3D 


<4-4>  st  ° 3V 


=0  (v  = inward  normal) 


permits  us  to  reduce  the  bilinear  form  H 


to  the  case  where 


only  the  boundary  data  of  the  first  kind 


_ + n-1  -t 

3 u 3 u 

— ^ = f , , ...  , = f on  3 D 

3 v 1 n-1  n 

3 v 


This  representation  is  more  easily  computed  using  the  complex 


notation.  To  this  end  we  recall  the  elementary  identities, 


J 
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A 1 

3k*mu 

(4.6) 

ds  1 

SzSz*1 

f <k  1 

1 d u 

k + 

1 k 
= 1 

j 

1 

idv\ 

I ^k+m+1 
(Tzk  + 13z  *ml  tiS 


Stl*  4S  . I -=~AX  -51 . 

3zk3z‘™  ds 


k 

l < 

l=o 


>k-‘(kJ 


,k 
3 u 


3zk_l3z*1 


k-2t 


which  hold  on  3D.  The  first  two  brackets  of  (4.2)  we  compute 
directly  as 


(4.7) 


where 


(4.8) 


(4.9) 


{Vt  Mlutl  - Ut  M* [Vt ] ) - (V  L ( U t ] - Ut  L* I V ] ) 


~ p + —5—  p* 

Dz  3z  * 


n k-1 


==  l l (vtAkV  k.t-i — T 

k=  1 £=o  3z  1 K3  3zK  *•  i3z*:) 


3k+3-t'1U 


l=o 
n-1  k-1 


, t 3 


- I l (-I)"  — ^VBki) 

k=l  t=o  3Z  3z 

j=o 


3k+j-l-l 


k_l"18z*j  ' 


- n . j-1  .m+k 

r*  :=  J'-11  J0('1)ra^sp?  ‘Vkj’ 

k~o  m=o  3 2 * 3z  J 

3 = 1 


g j-m-ljjt 

~ m 
3z 


n- 1 . j-1  m .m+k 

l (-D k i (-Dm  3 

k-o  m=o 


3J 


-m-  1 , 


a 4~m~k  lVBkj> 
3z*  3z 


3z 


1 = 1 

Then  (4.1)  takes  on  the  form 


(4.10) 
where 

(4.11) 


U ( z , z* , t ) = | dt 
o 


N (U,S)ds 


3D 


No(U'V)  :=  1 f 


15.  \ . 

ds  j 


This  in  turn  may  be  expressed  in  terms  of  tangential  and  normal 
derivatives  of  the  data  using  (4.6).  When  S is  a Green’s  function, 
then  (4.10)  is  directly  evaluated  in  terms  of  (4.5). 
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